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' In this paper we emphasize the importance of the Sloan Digital Sky Survey 

^ ■ (SDSS) QSO clustering statistics as a unique probe of the Universe. Because the 

^ . complete SDSS QSO sample covers a quarter of the observable universe, cosmological 
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parameters estimated from the clustering statistics have an implication as a test of 
the cosmological principle, by comparing with those from local galaxies and other 



00 I cosmological observations. Using an analytic approach to the power spectrum for the 



QSO sample, we assess the accuracy with which the cosmological parameters can be 



PsJ ■ determined. Arguments based on the Fisher matrix approach demonstrate that the 

^ ' SDSS QSO sample might have a potential to provide useful constraints on the density 

■ parameters as well as the cosmic equation of state. 
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1. Introduction 



The cosmological principle, the assumption of the homogeneity and isotropy of the Universe 
on large scales, is one of the most fundamental in the framework of the cosmology (e.g., 
Peebles 1993). As Lahav (2001,2002) reviewed, observations of the comic microwave background 
(CMB), the cosmic X-ray background, radio sources, and the Lyman-a forest strongly support 
the cosmological principle, though it is difficult to prove the principle definitely. Actually 
inhomogeneous cosmological models, which challenge the cosmological principle, have been 
proposed (Kantowski & Thomas 1998, Kantowski 2001, Barrett & Clarkson 2000, Celeria 2000). 
For example, it is claimed that a large local void universe might be viable, which is compatible 
with the observations of type la supernovae, the cosmic microwave isotropy, and the distribution 
of the local galaxies (Tomita 2000, 2001). 

We emphasize the importance of the QSO clustering statistics as a unique and independent 
probe of the Universe. Especially it provides a chance to test the cosmological principle by 
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comparing the cosmological parameters measured from the QSO clustering analysis with those 
from local galaxies and other cosmological observations. Recently the two degree field (2dF) QSO 
redshift (2QZ) Survey reported that a simply biased mass distribution explains the QSO spatial 
clustering and that the cosmological parameters can be measured (Groom et al. 2001, Hoyle et al. 
2002a; 2002b, Outram et al. 2001). Their preliminary results with the 10k catalogue favor the 
cold dark matter (CDM) model with a cosmological constant, though the constraint is not very 
tight (see also Yamamoto 2002; Hereafter Paper I). 

On the other hand, the Sloan Digital Sky Survey (SDSS) will aim to compile a homogeneous 
catalogue of 10^ QSOs as well as 10^ galaxies. The SDSS stores the 5 bands CCD imaging of 10^ 
deg^ in the North Galactic Gap and of three 2.5 x 90 deg^ strips in the South Galactic Gap. The 
QSO spectroscopic survey of the North Galactic Gap will detect 10^ QSOs in the range of the 
redshift ^ z ^ 5 with the limiting magnitude m = 19 ~ 20. A QSO catalogue with the SDSS 
early data release has been reported (Schneider et al. 2002). Because the SDSS QSO sample 
will cover a quarter of the observable Universe, the cosmological parameters measured with the 
clustering statistics provide a unique implication for the cosmological principle. Therefore it is 
important to estimate the accuracy with which the QSO clustering statistics can measure the 
cosmological parameters. 

Several authors have proposed possible methods (e.g., Ballinger, Peacock, & Heavens 1996, 
Matsubara &: Suto 1996, Popowski et al. 1998), which originate from the Alcock-Paczynski's 
test (1979). Very recently, Galvao et al. (2002) have claimed that a significant constraint on 
the density parameters as well as the cosmic equation of state of the dark energy, which is a 
quite important issue in physics and cosmology (see e.g., Caldwell, Dave &; Steinhardt 1998), can 
be obtained from the Alcock-Paczynski test with the 2dF QSO sample. Concerning the SDSS 
QSO sample, however, a systematic assessment on the capability of measuring the cosmological 
parameters have not been well investigated, as far as we know. The reason might come from the 
fact that physical processes of the QSO formation and time-evolution have not been completely 
understood. On the other hand, the number count of the SDSS QSO sample in the early data 
release has been demonstrated by Schneider et al. (2002), which we use in our investigation. 

Motivated by these situations, we revisit the QSO clustering statistics and assess the accuracy 
with which the cosmological parameters can be measured by the SDSS QSO sample.^ For this 
purpose we adopt the Fisher-matrix approach, which has been often used to assess the accuracy 
of measurement of the cosmological parameters from large surveys such as GMB anisotropy 
observations, galaxy redshift survey, and supernova data sets. The Fisher matrix approach bases 
the method on mathematical theorems, which enables us to clearly interpret results. The method 
adopted in the present paper is a simple extension of an approximate Fisher matrix approach 
using the power spectrum, which was originally formulated for galaxy samples (Tegmark et al. 
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1998; Tegmark 1997). We use the power spectrum analysis extended so as to incorporate the 

various observational effects, the hght-conc effect, the hncar redshift distortion and the geometric 
distortion, which is useful for QSO samples (Paper I, Suto ct al. 2000, Yamamoto & Suto 1999, 
Yamamoto &; Nishioka 2001). This paper is organized as follows: In section 2, we briefly review 
formulas to evaluate the Fisher matrix element. In section 3, the best statistical errors of the 
cosmological parameters are shown assuming the complete SDSS QSO sample, which demonstrates 
the usefulness of the power spectrum analysis of the sample. In section 4, the problem of 
degeneracy in the parameter space is briefly discussed. Section 5 is devoted to conclusions. 
Throughout this paper we use the unit in which the light velocity equals 1. 



2. Method 

The Fisher matrix approach provides a useful technique to estimate errors in measuring 
cosmological parameters from a given data set (see e.g., Tegmark, Taylor & Heavens 1997, 
Jungman et al. 1996;1997). Tegmark et al. (1998) developed a useful approximate method to 
estimate the Fisher matrix element using the power spectrum of large surveys of galaxies. The 
method used in the present paper is a simple extension of their work. In an analysis of the QSO 
clustering, the additional observational effects, the light-cone effect and the geometric distortion 
effect, must be taken into account. The power spectrum analysis incorporating these effects can 
be developed by extending the work by Feldman, Kaiser & Peacock (1994; hereafter FKP). We 
start from briefly reviewing this power spectrum analysis useful for the QSO samples. 

First we assume that the number density field of sources n(s) is constructed from a data 
set, where s denotes the three dimensional position of the sources. Here we assume that a 
distance-redshift relation s(= |s|) = s{z) was adopted to plot the map (see below for the explicit 
definition s{z)). Following FKP, we define the density fiuctuation field by 

_ n(s)-an,(s) 
[/dsn(s)2]V2' 

where a is a constant, ns(s) is a synthetic catalogue, n{s) is the expected mean number density. 
The optimum weight factor introduced in FKP can be fixed as a constant because nP{k) is less 
than unity for the QSO sample. This means that the power spectrum analysis considered in the 
present paper does not require information on the evolution of the power spectrum as a function 
of the redshift. The power spectrum is constructed by the square of the Fourier coefficient of the 
fluctuation field defined by 

J^(k) = j dsF(s)e*^■^ (2) 

Then the estimator of the power spectrum (before averaging over a shell of the radius k{= |k|) in 
the Fourier space) is defined by subtracting the shot noise contribution Pghot from |.^(k)| 

P{k) = iJ^ik)^ - P.hot. (3) 
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Using the two-point functions of n(s) and ns(s) (see equation (2.1.5) in FKP), we can compute 
the expectation value of (|.?^(k)p), which yields 

Jdz{dN/dz) 
J dzW{z 

and 



w = (! + «) -' , (4) 



_ JdzWiz)P^ik,z) 
JdzW{z) ' 

where we defined W{z) = {dN / dzY {dz / ds) / by using the observable quantity dN/dz, the 
number of sources per unit redshift and unit solid angle, instead of n(s), which are related by 
dN = dss^n(s), P'*(k, z) is the local power spectrum defined on a constant time surface of the 
redshift z. In deriving (5), we adopted the distant observer approximation, |si — S2I << |si + S2I. 
Under this approximation, the correlation function is related with P"(k, z) as 

^(si, S2) = (27r)-3 J dkP-{k, z{\si + S2|/2))e^''-(^i-^2). (6) 

In a similar way, the variance of the power spectrum (covariance matrix) can be approximately 
estimated 

(AP(k)AP(k')) = ^^^(5(k - k'), (7) 

K 

where AP(k) = P(k) - (P(k)), and 

K = AnJ dss^n{sf = AnJ dzW{z) (8) 

with the solid angle of a survey area AO. Equation (5) is the power spectrum including the 

light-cone effect, which is incorporated by averaging the local power spectrum P"(k, z) over the 
redshift with the weight factor W{z). The usefulness of the power spectrum (5) is discussed by 
comparing with results of the 2QZ group (Paper I). 

In general, the Fisher matrix element is defined by 

-d^ InL 



~ \ de.dOj /' 

where 9i denotes the parameters and L is the probability distribution function of a data, given 
the model parameters. For simplicity, we adopt the approximation of the Gaussian probability 
distribution function for AP(k) 

L oc exp[-i J dkj dk'AP(k)C(k,k')"^AP(k')], (10) 

where the matrix C(k, k')"-*^ is the inverse of the covariance matrix (AP(k)AP(k')). Following 
the work by Tegmark (1997), the Fisher matrix element reduces to 

F Y (21+1) r-^'^dkk^^^^^^^ 



=0,2, 



(11) 
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where 

{Pi{k))= [' dfiCi{i,){P{k)), (12) 
Jo 

and is the Legendre polynomial of the Z-th order, and /x denotes the directional cosine 

between the line-of-sight direction and the wave number vector. Note that (-P(k)) is a function 
of k and fi under the distant observer approximation. This allows P"(k, 2) in equation (5) to be 
written P"'(k, /i, z), which we model as (Ballinger et al. 1996) 

P«(A;,/x,s) = -\ 

C||C^ 

^PQso[q\\~> — ,gx^ ,z), (13) 

where Pqso{Q\\iQ±i z) is the QSO power spectrum, and are the wave number components 
parallel and perpendicular to the line-of-sight direction associated with the comoving coordinate of 
the real universe, Cy and Cj_ are defined by Cy = dr{z)/ds{z) and Cj_ = fK{r{z))/s{z), respectively, 
with the comoving distance r{z) and the angular diameter distance fxiriz)) of the real universe. 
The explicit expression of r{z) is presented by equation (15). We model the QSO power spectrum 
of the spatial distribution incorporating the linear distortion (Kaiser 1987, see also Yamamoto 
Nishioka & Suto 1999), as follows: 

PQSo{qii,q±,z) = (^l + ^^yb{z)''P^^,{q,z), (14) 

where f{z) = d\nD{z)/d\na{z) with the linear growth rate D{z) and the scale factor a{z), 

= qfl + qIj Kz) is the scale independent bias factor, and PmassiQ, z) is the CDM mass power 
spectrum. The term in proportion to f{z) describes the linear distortion due to the peculiar 
velocity on large scales. We work within the linear theory of density perturbations because 
we consider the large scale QSO clustering statistics. In the present paper, we assume the 
Harrison-Zeldovich initial mass power spectrum and adopt the fitting formula of the cold dark 
matter transfer function by Eisenstein and Hu (1998), which is robust even when the baryon 
fraction is large. 

Equation (11) is an extension of the previous result by Tegmark (1997). In the limit that the 
linear distortion, the geometric distortion, and the light-cone effect are switched off, equation (11) 
reproduces the result by Tegmark (1997). In the investigation in the present paper, the light-cone 
effect and the geometric distortion effect are particularly important. For a sample of a narrow 
range of the redshift, the light-cone effect can be negligible, however, it should be properly taken 
into account in the QSO spatial power spectrum analysis of the 2QZ sample and the SDSS 
sample, which consist of sources in the wide range of the redshift. In general, the geometric 
distortion becomes substantial for high-redshift samples. On the other hand, the linear distortion 
effect is rather minor in our investigation. It is well recognized that the linear distortion effect 
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is substantial and problematic in the Alcock & Paczynski's geometric test (e.g., Ballinger et al. 
1996). It is the fact that the anisotropic component of the power spectrum is sensitive to the 
linear distortion effect even at the high rcdshift. However, as we show in the next section, the 
dominant contribution constraining the cosmological parameters comes from the isotropic part 
of the power spectrum. The linear distortion effect only increases the overall amplitude of the 
isotropic part of the power spectrum in the scale independent bias model. This is the reason why 
the linear distortion is the minor effect in our investigation. This situation is in contrast to case 
of low redshift galaxy samples, in which the light-cone effect and the geometric distortion are 
negligible, while the linear distortion is rather important. 



3. Results 

In general, the comoving distance is written, 

''''^^^Jo Ho[n„,{l + z'f + QKil + z'y + 0,(1 + z')3(i+-q)]1/2 ' 

where the Hubble parameter is Hq = 100 /ikm/s/Mpc, ^l„i, 17/^, and Qx arc the cosmological 
parameters of the densities and curvature, two of which are independent by virtue of the relation 
Clra + + = ^1 and we introduced the constant cosmic equation of state parameter, wq. In 
the present work, we focus on the four independent parameters ilm, ^^6, and wq. Here ilm is 
the sum of the CDM and the baryon components. Following an analysis by the 2QZ group (Hoyle 
et al. 2002a), we choose s{z) to be the comoving distance with $7^ = 0.3, CIk = 0, and wq = —1: 

^ Jo i?o[0.3(l + 2')' + 0.7]V2- 

Our results does not significantly depend on this choice of s{z). For dN/dz in the present paper, 
we adopt the recent result reported by Schneider et al. (2002) with the first SDSS QSO catalog. 
The QSO clustering bias in the SDSS sample has not been well investigated at present. However, 
Groom et al. reported on the time-evolution of the QSO bias by analyzing the 10k catalog of the 
2QZ Survey (2001). According to them, the clustering amplitude does not show significant time 
evolution. For modeling the bias, following (Calvao et al. 2002), we first adopt the form 

h{z) = 1 + (6o - l)D{z)-'-\ (17) 

and determined the constant 6o(= 1-2) so that the power spectrum best matches the 2QZ 
result (Hoyle et al. 2002a). For the linear growth rate of the model with wq ^ —1, we adopted 
the fitting formula developed by Wang and Steinhardt (1998). Thus we vary the linear growth 
rate depending on cosmology. However, our result is not significantly altered by this modelling 
because the linear distortion effect is not substantial as described in the previous section. 

1/2 

Figure 1 shows the best statistical errors, A^j = 1/ -^j/ j some of which are normalized by 
the target parameters (the parameters supposed to be true), = 0.28, $7^ = 0, Jib = 0.04 and 
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WQ = —1, as functions of kmax in equation (11). We fixed kmm = 0.005 hMpc and assumed tt 

steradian of the survey area. The sohd and dotted curves in each panel arc the I = and I = 2 
components of Fij, respectively. The dashed curve corresponds to the case that both components 
are summed up. Contribution of the higher modes I > 2 to the Fisher matrix element is small. 
Thus the contribution constraining the cosmological parameters is dominated by the component 
I = 0. The component I = 2 corresponds to the constraint from the anisotropic power spectrum 
{P2{k)), while I = corresponds to the constraint from the isotropic power spectrum (Po(fc))- The 
geometric distortion effect causes anisotropics of the power spectrum, and the Alcock-Paczynski's 
geometric test utilizes this fact. The geometric distortion effect also alters the shape of the power 
spectrum {Po{k)) due to the scaling of the wave number, equation (13). Therefore the constraint 
on the cosmological parameters in the present paper mostly relies on the latter effect that the 
shape of the isotropic power spectrum changes due to the geometric distortion. 

1/2 

Figure 2 shows the best statistical errors, A9i = l/F--' , as function of /cmin with fixed 
^max = 1 hMpc^^. The target parameters and the meaning of the curves are same as those 
of Figure 1. These two figures indicate that the Fourier modes of k ^ O.Ol/iMpc"^ and 
k ^ 0.1/iMpc~^ are not significantly relevant to the integration for the Fisher matrix elements. 
Thus the Fourier modes of 0.01/iMpc^^ ^ k ^ 0.1/iMpc^^ dominate over the constraint on the 
cosmological parameters. This is a contrast to the case of the galaxy sample (Tegmark 1997). This 
difference traces back to the shot-noise dominance in the QSO sample. This result also suggests 
the validity of our assumption that nonlinear effect of density perturbations is small because the 
nonlinearity is influential on scales k ^ O.lhMpc"^ . 

In the present paper, we have not considered the effect of the window function depending on 
the shape of a survey volume. The effect of the window function on the power spectrum for the 
2QZ sample was investigated in the literature (Hoyle et al. 2002a). According to their result, 
the 2QZ power spectrum might be affected for k ^ 0.02/iMpc~^. Though a definite investigation 
using a mock sample is needed, for the case of the SDSS QSO sample, it might be expected that 
the eflFect of the window function appears on the larger length scales (smaller k) than that of the 
2QZ power spectrum. From figures 1 and 2, as long as the window function effect in the SDSS 
QSO sample appears on the scale k ^ O.Ol/iMpc"^, the estimation of the errors is not altered 
significantly. 

We note that our result might depend on the modelling of the bias. This is because the Fisher 
matrix element is in proportion to the amplitude of the power spectrum of the QSO distribution. 
Figure 3 shows the same as Figure 1, but with adopting the bias model b{z) = bQD{z)^^, in which 
the amplitude of the QSO clustering is almost constant as a function of the redshift z. Here we 
set bo = 0.8 so as to best match the 2QZ power spectrum (see also Paper I). The best statistical 
errors are worse than those of Figure 1, though the difference is not large. The reason of the 
smallness of the difference is explained as follows: We normalized the power spectrum so as to 
match the 2QZ power spectrum, in which the 2dF QSO sample is distributed in the range of 
redshift 0.3 < z < 2.3. The difference is the contribution from the higher redshift z ^ 2. The 
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bias model (17) has larger amplitude at the higher redshift, which yields larger amplitudes of the 
power spectrum and the Fisher matrix elements. 

In the above investigations, we considered the simple bias models, whose amplitude we 

determined so that the power spectrum best matches the 2QZ result. This method might 
estimate the errors smaller. Then wc also investigated the best statistical errors marginalized over 
parameters of the bias. We here consider the bias model parameterized by bo and v, as follows: 

= l + (18) 

then marginalize the probability function with respect to bo and v. Explicitly, we considered the 
space of the parameters {Qm, bo, v) or {wq, bo, u) , and consider the Fisher matrix element 
marginalized by integrating the probability function with respect to 6o and v. This process 
actually reduces the accuracy with which the cosmological parameters can be determined. For 
example, in the case with the target parameters 6o = 1-2 and u = 1.7, we obtained the relative 
error AJ7„,/f7,„ 0.4 and Awq ~ 0.7, where we assumed the same target parameters for the 
cosmological models as those of Figure 1. Similarly, we have Ar2„j/r2m — 0-35 and Aiuq ~ 0.5 for 
the target parameters bo = 1.5 and v = 1. We also have AO^/^m — 0.25 and Awq ~ 0.35 for 
6o = 1-5 and u = 1.7. 



4. Degeneracy 

In general the feasibility of determining the cosmological parameters is restricted by 
degeneracies in the parameter space. Actually many possible parameters can appear in generic 
cosmological model. Here we demonstrate some aspects of the degeneracies in the cosmological 
parameters. In Figure 4 we plot two dimensional confidence contours, where wc assumed the same 
QSO sample and the bias model as those of Figure 1, Vim = 0.28, IIa = 0.72, = 0.04, wq = — 1, 
and the Harrison-Zeldovich spectrum n = 1. Each panel shows the confidence contours on the 
two parameters, i.e., and the other parameter labeled in each panel, with assuming the rest 
parameters to be fixed. Thus this figure shows the degeneracy between Jl^ and other parameters. 
From this figure it is clear that the uncertainty about fi^ is problematic when determining 
It also shows that the degeneracy between $7^ and n, the index of the initial power spectrum, 
is relatively large. Thus the degeneracy between these parameters will be problematic in terms 
of determining all the parameters simultaneously. Nevertheless a useful constraint will be 
obtained from the QSO sample by combining it with constraints from Big Bang Nucleosynthesis, 
CMB anisotropics, la supernova, weak gravitational lensing and other observations. Then we 
can emphasize that the QSO sample provide a chance of a worthy complimentary test for the 
cosmological models. 
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5. Discussions and Conclusions 

In this paper we have assessed the feasibihty of measuring the cosmological parameters 
by means of the power spectrum analysis of the SDSS QSO sample, with an emphasis on a 
unique observational probe of the Universe. If tight constraints are provided, the result will have 
implications for testing the cosmological principle and inhomogeneous cosmological models. The 
best statistical errors are estimated using the Fisher matrix approach, which is a simple extension 
of the previous works for galaxy samples. For example, the best statistical errors are at 10% level 
for fljn and Qk, a few times 10% for wq and O^. These errors can increase to several xlO%, 
depending on the bias amplitude, when marginalizing over the parameters of the bias. Thus our 
result shows that the SDSS QSO sample might have a useful potential as a test of cosmological 
model. But, this does not mean that we have shown the capability of the Alcock-Paczynski test 
using the anisotropics of the power spectrum, because the constraint is almost from the isotropic 
component of the power spectrum. The cosmic degeneracy limits the capability of the method 
in determining all parameters with the best statistical accuracy simultaneously. However, the 
point is the uniqueness of the method as a probe, and the result would be useful for more robust 
understanding of the Universe. 

Finally we mention comparisons with other works related with the present paper. Ballinger 
et al. (1996) used a method similar to the present paper, to predict a constraint from the 2dF 
QSO clustering statistics. Outram et al. (2001) have applied their method to the 2QZ sample. 
The 1-cj error of Q„i by Outram et al. (2001) is consistent with our result obtained by applying 
our method to the 2QZ sample. Very recently, Matsubara and Szalay (2002) have also reported 
a similar investigation based on an eigenmode analysis, while our method relies on the power 
spectrum analysis. Concerning the QSO sample, as we have discussed, the number density in 
proportion to dN/dz and the modelling of the bias b{z) can be sensitive to the final results. These 
modellings and the choice of the cosmological redshift space s{z) are also different. 
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Fig. 1. — The best statistical errors A9i{= 
normahzed by the target parameters = 0.28, = 0, = 0.04 and wq = — 1. The sohd and 
the dotted curve in each panel corresponds to I = and / = 2 component, respectively. The dashed 
curve is the sum of the both. We adopted the solid angle of the survey area, AO, = tt steradian. 
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Fig. 2. — The best statistical errors A9i{= 1/F-- ) as function of fcmin- The target parameters and 
the meanings of the curves are same as those of Figure 1. Here we fixed femax = 1 /iMpc~^. 
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Fig. 3. — Same as Fig. 1 but with the bias model of the form b{z) = boD{z) ^, where 6o(= 0.8) is 
used. 
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Fig. 4. — Confidence regions. We have chosen the target parameters as Qrn = 0.28, U\ = 0.72, 
rifo = 0.04, WQ = — 1 and n = 1. In the panels, we used the superscript 'tr' to denote the target 
parameters explicitly. Each panel shows the confidence contours on the two parameters, and 
each of the other parameters, which is labeled in the panel, assuming that the rest parameters 
are fixed. The curves are the contours of the 68%, 95% and 99.7% confidence contours. Axes are 
normalized by the target parameter. 



